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SIMULATION  OF  COMPLEX  SHOCK  REFLECTIONS 
FROM  WEDGES 'IN  INERT  AND  REACTIVE  GASEOUS  MIXTURES 


1.  Advances  in  FCT  Techniques 

In  this  paper  we  describe  new  adaptations  of  the  Flux-Corrected  Transport  (FCT) 
algorithms  developed  by  Boris  and  Book  (1973,  1976)  for  solving  fluid  equations,  and 
discuss  their  application  to  multidimensional  shock  reflections  in  inert  and  reacting 
gaseous'mixtures.  In  particular,  we  consider  planar  constant-velocity  shocks  reflecting 
from  wedges.  Under  certain  circumstances  double  Mach  stems  are  formed.  Historically 
these  have  proven  to  be  difficult  to  calculate  with  high  accuracy,  although  many  schemes 
have  been  available  to  analyse  compressible  flow  on  a  computationally  discretized  mesh: 
the  method  of  characteristics,  spline  techniques,  Glimn-type  random  choice  schemes,  and 
finite  element,  finite  difference,  and  spectral  methods.  We  believe  that  the  calcula- 
t tonal  difficulties  experienced  on  this  problem  were  the  result  of  excessive  numerical 
diffusion,  especially  in  the  region  of  the  contact  surface. 

In  any  Ruler lan  calculation,  numerical  diffusion  arises  because  material  which  has 
Just  entered  a  computational  cell,  and  is  still  near  one  boundary,  becomes  smeared  over 
thb  whole  cell.  FCT  minimises  this  effect.  FCT  algorithms  can  be  constructed  as  a 
weighted  average  of  a  low-order  and  a  high-order  finite-difference  scheme.  If  the  fluid 
equations  are  written  in  conservative  form,  both  schemes  are  Implemented  using  trans¬ 
port  ive  fluxes.  Rach  flux  describes  the  transfer  of  mass  for  some  other  extensive 
quantity)  from  one  point  to  a  neighboring  point.  The  procedure  for  assigning  weights 
involves  limiting  or  "correcting"  the  fluxes  at  certain  points.  The  higher-order  scheme 
is  used  to  the  greatest  possible  extent,  consistent  with  .welding  the  introduction  of 
dispersive  ripples  (undershoots  and  overshoots).  The  weights  for  the  low-order  scheme 
are  chosen  to  be  Just  sufficient  to  eliminate  these  ripples,  thus  assuring  the  property 
of  "monotoniciiy"  or  "positivity."  The  result  is  an  algorithm  which  effectively  reduces 
to  the  higher-order  scheme  wherever  the  fluid  properties  change  gradually.  Sear  sharp 
discontinuities,  however,  enough  diffusion  is  supplied  to  retain  nenotpnleity.  At 
shock  fronts  this  procedure  automatically  produces  the  correct  local  viscous  heating. 


The  prototype  second-order  finite-difference  formula 
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illustrates  the  procedure.  Mere  j  labels  grid  position,  n  denotes  time  level  and 
t.+,  *v.  .At/#*.,,  and  V. are  dimensionless  advert  ion  and  diffusion  coefficients, 
respectively.  JWd  vriteJ*N  S  |+c),  where  e  is  a  "dipping  factor"  measuring  the 

extra  diffusion  added  to  achieve  positivity.  When  e  •  0,  the  above  scheme  is  second- 
order;  in  the  vicinity  of  shocks  c  '**.  1  and  it  effectively  reduces  to  first  order. 


A  numerical  diffusion  Reynolds  number  Ole) 
characterist  ic  site  of  a  structure  in  tne 


2L/cAs  can  be  defined,  where  t»  is  the 
Kven  the  most  accurate  spectral  simu¬ 


lations  require  setting  c  »  1  to  guarantee  positivity  linearly.  This  gives  rise  to  the 


usual  definition  of  the  numerical  Reynolds  number,  21./ A*.  Algorithms  such  as  FCT  which 
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guarantee  monotonicity  nonlinearly  can  have  average  values  \c'  'v  10  — 10  ,  intro¬ 

ducing  much  less  overall  dissipation  and  permitting  calculations  with  effective  Reynolds 

numbers  such  that  Re  "v(Re),m  >>2L/6x. 

ND 

Four  advances  in  FCT  techniques  have  enabled  us  to  perform  a  series  of  shock  and 
detonation  calculations  with  high  accuracy.  These  techniques  are  easy  to  program,  and 
they  have  wide  applicability  to  general  quasi-linear  hyperbolic  equations  (i.e.,  equa¬ 
tions  describing  continuum  conservation  laws). 

The  first  of  these,  a  generalization  of  FCT  due  to  Zalesak  (1979),  removes  the 
necessity  of  tinestep  splitting  in  multidimensional  hydrodynamics.  This  reduces  errors 
associated  with  time  splitting  in  regions  of  the  flow  which  are  nearly  incompressible. 

The  second  refers  to  the  development  of  FCT  algorithms  in  which  the  spatial  derivatives 
can  be  approximated  to  arbitrarily  high  order  (fourth,  sixth,  eighth,  etc.,  or  pseudo- 
spectral).  These  innovations,  which  relate  to  the  transport  algorithm  itself,  have 
been  implemented  in  a  two-dimensional  hydrocode  which  utilizes  the  leafrog-trapezoidal 
(L-T)  algorithm  and  is  therefore  dissipationless.  Both  complex  and  double  Mach  stem 
structures  are  obtained  (cf.  Ben-Dor,  1978}  Ben-Dor  and  Glass,  1978,  1979), 

The  third  new  technique,  adaptive  rezoning,  is  an  extension  to  two  dimensions  of 
the  dynamic  rezoning  employed  in  detailed  one-dimensional  reactive  flow  simulations  by 
Oran,  ec  al,  (1979).  This  concentrates  needed  spatial  resolution  in  the  vicinity  of 
moving  shocks,  contact  discontinuities  and  reactive  surfaces.  The  technique  is  illus¬ 
trated  with  shock  calculations  using  a  time-split  code  (FAST2D).  In  air  for  M=5  and 
6  ^S8,  the  results  fall  very  close  to  the  boundary  between  regular  and  Mach  reflection. 
The  calculated  wall  pressures  are  in  detailed  agreement  with  the  results  of  Bertrand 
(1972).  The  fourth  technique  is  a  generalization  of  the  induction  time  approximation 
used  in  earlier  flame,  ignition  and  shoek  work  (Oran  et  al,  1980a,  b).  This  provides  a 
simple,  efficient,  yet  reasonably  accurate  global  chemical  kinetics  package  to  be  used 
in  connection  with  these  comprehensive  two-dimensional  hydrodynamics  calculations. 

Section  2  describes  the  results  of  calculations  in  which  a  planar  shock  is  reflected 
from  a  wedge  in  an  inert  gas.  In  Section  3  we  present  the  results  of  calculations  of 
detonations  initiated  by  shock  reflections  in  stoichiometric  mixtures  of  II,  in  air  at 
low  pressure.  Section  4  summarises  our  conclusions.  2 

2.  Shock  Reflections  in  Air 

The  utility  of  these  advanced  FCT  methods  has  been  demonstrated  by  applying  them 
to  transient  reflections  of  planar  shocks  from  wedges  for  various  shock  strengths  M  and 
wedge  angles  8  •  For  nonreaeting  flews  at  Mach  numbers  greater  than  about  2.3  and  wedge 
angles  between  20  and  50  degrees,  double  Mach  stems  can  develop.  Numerical  schemes 
previously  used  for  this  problem  reproduce  qualitatively  the  wave  structure  and  shape, 
but  have  difficulty  making  accurate  predictions  of  flow  details  such  as  density  contours 
(a  conclusion  drawn  by  Bea-ttar  and  Glass,  1978)  even  in  the  single  Mach  stew  case.  To 
our  knowledge,  successful  calculations  of  the  double  Mach  stem  ease  have  not  vet  been 
published.  In  this  paper  we  discuss  the  series  of  calculations  summarized  in  Table  1. 

Open  boundary  conditions  are  used  on  the  left,  right,  and  top  edges  of  the  mesh, 
i.e.,  density,  pressure  and  velocity  are  set  equal  to  their  pre-shock  or  post-shock 
values,  depending  on  whether  the  incident  shock  front  has  passed  that  point.  Reflecting 
conditions  are  imposed  on  the  bottom  of  the  mesh,  which  corresponds  to  the  wedge  surface. 
Examples  of  the  calculated  density  contours  and  wave  structure  for  the  double  and  complex 
Mach  reflection  cases  are  shown  in  Fig.  1.  The  Incident  shock,  I,  the  contact  surface, 
CS,  and  the  first  and  second  Mach  stems,  M,  and  )L,  are  indicated  in  Fig.  la,  Note  In 
particular  the  forward  curl  of  the  contact ‘surface  near  the  wail  and  the  small  region 
(4  by  7  mesh  points)  of  high-density  gas  just  to  the  left  of  tin  point  where  the  contact 
surface  impacts  the  wall.  The  latter  causes  a  second  peak  *  .the  pressure  and  density 
distribution  on  the  wall,  as  shown  in  Fig.  lc.  The  accuracy  a.  Jut astuiatlons  has  been 
verified  by  comparison  with  experimental  density  distributions  along  the  wall,  as  shown 
in  Fig.  2,  and  with  experimental  pressure  measurements  (Bertrand,  1972).  Note  that  FCT 

Srovidex  adequate  resolution  of  the  key  surfaces  (Contact  surface  and  second  Mach  stem) 
n  regions  as  small  as  S  by  S  cells. 
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Two  additional  cases  were  calculated  with  larger  values  of  0  (Fig.  3).  As  the 
wedge  angle  increases,  the  tlach  stem  develops  fiiore  slowly,  being  separated  from  the 
wedge  by  a  triple-point  angle  of  only  one  or  two  degrees.  [The  triple-point  angle  a 
is  the  angle  subtended  by  the  Mach  stem  as  viewed  from  the  end  of  the  wedge  at  which 
the  shock  was  first  incident  (Fig.  4.)  To  reduce  the  size  of  the  mesh  needed,  it  is 
necessary  to  calculate  in  the  frame  of  the  Mach  stem  and  to  rezone.  For  these  calcula¬ 
tions  we  employed  the  time-split  code  FAST2D  (with  a  130  x  50  mesh),  which  incorporates 
an  automatic  continuous  ("sliding  zone")  regridding  procedure  (Oran  et  al.,  1979).  For 
the  small  0  cases  discussed  above,  where  regridding  is  necessary,  FAST2D  yielded  re¬ 
sults  very  similar  to  those  obtained  with  the  L-T  code.  The  cases  with  wedge  angles 
of  44°  and  46.5°  constitute  a  severe  test  of  the  numerical  algorithm  because  of  the 
small  triple-point  angle.  Because  a  is  approximately  equal  to  2.8°  and  1.5°,  respec¬ 
tively,  considerable  spatial  resolution  and  a  large  amount  of  running  time  are  usually 
necessary  to  get  accurate  flow  fields.  Figure  4  illustrates  our  adaptive  rezone  tech¬ 
nique  on  a  grid  of  60  x  40  cells  with  varying  cell  dimensions,  5x,  6y.  This  method 
requires  one-fifth  the  number  of  cells  required  in  a  uniform  grid  calculation.  A 
uniform  region  consisting  of  the  smallest  cells  covers  part  of  the  incident  shock 
front,  the  Mach  stems,  and  the  reflected  shock  structure.  Outside  this  finely  gridded 
region,  we  have  transition  zones  in  which  the  cell  dimensions  increase  smoothly  to  their 


maximum  values,  and  106y 


'min’ 


We  have  investigated  the  accuracy  of  the  numerical  simulation  by  comparing  the 
results  with  experimental  data  (Bertrand  et  al,  1972).  Because  the  eases  0  ■  44* 
and  46.5°  are  so  similar,  we  will  discuss  only  the  latter.  The  computed  value  of  a 
for  0  =  46.5®  is  approximately  2.5®  for  a  real-air  equation  of  state  and  approximately 
3.2®  tor  an  ideal  gas  with  y  ■  t.35.  Both  agree  with  the  measurements  to  within  the 
experimental  errors,  u  21.  In  Fig.  5,  we  compare  the  calculated  (using  a  real-air 
equation  of  state)  and  experimental  values  of  the  pressure  at  the  surface  of  the  wedge 
for  t1  «  46,5®.  The  agreement  in  the  shapes  of  the  pressure  curve  is  striklng—the 
values  of  the  lower  pressure  peak,  corresponding  to  the  Mach  stem,  are  nearly  identical. 
The  calculated  value  of  the  second  pressure  peak  is  11?  lower  than  the  experimental 
value  and  Is  thus  within  experimental  uncertainty.  Figure  6  shows  the  history  of  the 
pressure  on  the  wedge  calculated  for  an  ideal  gas  with  y  «■  1.35.  We  note  that  the  curves 
have  much  the  same  shape  as  for  the  real-air  simulations?  however,  the  first  pressure 
peak  is  again  117  lower  than  the  experimental  value.  Figure  J  shows  that  the  double 
Mach  reflection  shock  structure  is  well  resolved  in  the  simulation. 

3.  Detonations 

We  have  also  considered  analogous  shock  reflections  in  reactive  gases  (stoichio¬ 
metric  mixtures  of  IL  in  air)  at  low  pressure  (0.1  atm).  The  Induction  time  hypothesis 
(e.g.,  Oran  et  al,  3980)  represents  the  chemistry  through  a  composite  process,  in  which 
reactants  begin  to  combine  into  combustion  products  only  after  a  finite  time  has  elapsed. 
The  rate  at  which  the  energy-releasing  react  ions  proceed  depends  upon  a  single  parameter  , 
the  induction  time.  This  in  turn  is  a  function  of  the  local  thermodynamic  variables. 


Figure  8  shows  the  time  development  of  a  detonation  initiated  by  a  weak  reflecting 
shock.  The  incident  shock  was  chosen  so  that  the  pressure  behind  it  is  too  low  to  cause 
detonation  to  take  place  within  the  time  of  the  calculation.  As  with  the  calculation 
of  Section  2,  we  have  used  open  boundary  conditions  at  the  sides  and  top  of  the  system. 
The  sequence  of  six  pressure  contour  plots  traces  the  evolution  of  a  detonation  wave 
initiated  by  complex  Mach  reflection  at  the  surface  of  the  wedge.  Figure  8a  shows  the 
pressure  contours  corresponding  to  an  incident  shock  with  «  *  23*  and  M  *  4.0  which 
has  Just  begun  to  reflect.  The  Mach  stem  is  initially  too  smalt  to  be  resolved.  In 
Fig.  St»  the  Marti  stem  becomes  discernible,  but  as  yet  no  apparent  reaction  has  occurred. 
8y  frame  (c)  the  material  lias  begun  to  Ignite  at  a  position  well  behind  the  current 
location  of  the  Mach  stem.  When  the  Maclt  stem  passed  that  position  the  pressure 
increase  heated  the  mixture  sufficiently  to  cause  ignition  after  a  short  induction 
time  characteristic  of  the  IL'—air  mixture.  In  frames  (d)  and  (e),  at  later  times, 
the  pressure  at  the  Mach  stem  continues  to  grow,  leading  to  shorter  characteristic 
induction  times  for  material  between  the  Mach  stem  and  the  original  ignition  (mint. 

Thus  we  see  the  ignited  region  accelerate  along  the  wedge  surface  toward  the  .lech  stem, 
because  more  energy  Is  being  released  as  the  burning  continues,  the  boundary  of  the 
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ignited  region  also  accelerates  in  the  direction  of  the  reflected  shock  front,  com¬ 
pressing  and  heating  the  material  into  which  the  burned  gases  expand.  In  the  last 
frame  the  burn  front  has  overtaken  the  reflected  shock  and  Mach  stem,  as  we  see  from 
the  decrease  In  the  separation  of  the  pressure  contours  near  both  locations.  A  stable 
detonation  pattern  has  not  yet  emerged,  however.  This  is  evident  from  the  bending  out 
of  the  Mach  stem  and  the  lower  density  of  contours  between  the  Mach  stem  and  the  re¬ 
flected  shock/detonation  front. 

We  anticipate  that  shock  tube  experiments  will  confirm  this  wave  structure  and 
that  such  reactive  flow  calculations  will  be  extremely  useful  In  quantitatively  ex¬ 
plaining  the  experimentally  observed  multicell  structure  of  detonations  (Oppenheim  1970). 

4.  Conclusions 

Our  calculations  of  complex  and  double  Mach  reflection  are  in  close  agreement 
with  measurements  for  shocks  reflecting  in  air  from  wedges.  Because  of  the  accuracy 
and  speed  of  FCT  algorithms  and  the  effectiveness  of  adaptive  rozontng,  the  calculations 
are  accurate  and  economical  even  when  the  Mach  stem  develops  very  slowly.  All  of  the 
Important  features  (location  of  surfaces  of  discontinuity,  pressure  loading  on  the 
wedge  surface,  density  contours)  are  correctly  predicted.  The  results  do  not  depend 
sensitively  on  whether  the  l.-T  or  KAST2D  code  is  used.  Of  the  advances  discussed  in 
Section  1,  multidimensional  flux  limiting  and  the  adaptive  regridding  technique  seem 
to  he  the  most  efficacious  for  reflections  In  nonreactive  media.  We  conclude  that  FCT 
algorithms  reduce  numerical  diffusion  dramatically,  assuring  qualitative  improvements 
in  accuracy.  Uo  believe  that  to  achieve  comparable  accuracy  and  efficiency,  other 
hyd roe odes  must  employ  similar  nonlinear  algorithms  and  rc zoning  techniques. 

Our  calculations  in  reactive  gas  mixtures  show  that  detonations  tend  to  begin 
where  a  secondary  pressure  peak  arises  as  the  slip  surface  approaches  the  wedge. 

Because  of  the  finite  induction  time  in  our  kinetics  model,  the  detonation  begins 
somewhat  behind  this  pressure  peak.  The  high  resolution  our  calculations  achieve 
enables  us  to  follow  multiple  reflections  and  is  capable  of  providing  quantitative 
predictions  of  detonation  phenomena. 
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DOUBLE  MACH  REFLECTION 
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Fig,  1  -  (a)  Wave  structure  and  density  contours  for  double  Mach 
reflection  fro»  a  wedge;  (b)  reduced  density  contours  (  p/p )  for 
complex  Mach  reflection  with  levels  chosen  to  agree  with  tnose 
of  Ben-Dor  and  Glass  (1978);  (c)  corresponding  pressure  and  den¬ 
sity  profiles  on  the  wedge,  plotted  against  cell  number. 


6 


17 


ew  =  20° 

y  a  1.35 


Fig.  2  -  Comparison  of  the  calculated  density  profile  on  the 
wedge  (Fig.  1(h))  with  measure  *  values  (Ben -Dor  and  Glass,  1978 
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Fig.  4  -  Gridding  for  complex  shock  reflection  problems  with 
FAST2D  code.  Shown  are  incident,  reflected  and  Mach  shocks 
(solid  lines)  and  slip  surfaces  (dotted  lines)  for  incident 
shock  coming  from  the  left. 
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Fig.  5  -  Upper  and  lower  diagrams  show  pressure  in  PSIA  on  the 
wedge  as  a  function  of  position  for  o  -  46.5*(real  air  equa¬ 
tion  of  state)  at  two  times  in  the  simulation  and  as  a  function 
of  time  at  two  stations  In  the  experiment ,  respectively. 
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Albuquerque,  NH  87119 
Olcy  Attn  R  Hensley 

Boeing  Co. 

P  0  Box  3707 

Seattle,  WA  98124 

Olcy  Attn  H/S  42/37  R  Carlson 
Olcy  Attn  Aerospace  Library 

California  Research  &  Technology,  Inc. 

6269  Variel  Avenue 

Woodland  Hills,  CA  91364 
Olcy  Attn  Library 
Olcy  Attn  K  Krcyenhagen 

California  Research  &  Tech,  Inc. 

4049  First  Street 

Livermore,  CA  94550 
Olcy  Attn  D  Orphal 


General  Electric,  Co 
Space  Division 
Valley  Forge  Space  Center 
P  0  Box  8555 
Philadelphia,  PA  19101 
Olcy  Attn  H  Bortner 

General  Electric  Co. -Tempo 
816  State  Street  (P  O  Drawer  QQ) 
Santa  Barbara,  CA  93102 
Olcy  Attn  DASIAC 

Ceneral  Research  Corp. 

Santa  Barbara  Division 
P  0  Box  6770 
Santa  Barbara,  CA  93111 
Olcy  Attn  B  Alexander 

Higgins,  Auld  Association 
2601  Wyoming  Blvd  NE 
Albuquerque,  NH  87112 
Olcy  Attn  J  Bratton 
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IIT  Research  Institute 
10  W  35th  Street 
Chicago,  IL  60616 

Olcy  Attn  Documents  Library 
Olcy  Attn  R  Welch 
Olcy  Attn  M  Johnson 

Information  Science,  Inc. 

123  W  Padre  Street 
Santa  Barbara,  CA  93105 
Olcy  Attn  W  Dudziak 

Institute  for  Defense  Analyses 
400  Army-Navy  Drive 
Arlington,  VA  22202 

Olcy  Attn  Classified  Library 

J  H  Wiggins  Co.,  Inc. 

1650  S  Pacific  Coast  Highway 
Redondo  Beach,  CA  90277 
Olcy  Attn  J  Collins 

Kaman  Avidyne 
83  Second  Street 
Northwest  Industrial  Park 
Burlington,  MA  01803 
Olcy  Attn  Library 
Olcy  Attn  E  Criscione 
Olcy  Attn  N  Hobbs 
Olcy  Attn  R  Ruetenik 

Kaman  Sciences  Corp. 

P  O  Box  7463 

Colorado  Springs,  CO  80933 
Olcy  Attn  F  Shelton 
Olcy  Attn  Library 

Kaman  Sciences  Corp. 

Southern  California  Operations 
101  Continental  Blvd  Suite  85S 
El  Scgundo,  CA  9024 5 
Olcy  Attn  D  Sachs 


Lockheed  Missiles  &  Space  Co.,  Inc. 

P  0  Box  504 
Sunnyvale,  CA  94086 

Olcy  Attn  TIC-Library 

Los  Alamos  National  Scientific  Lab. 

Mail  Station  5000 

P  0  Box  1663 

Los  Alamos,  NM  87545 

Olcy  Attn  MS  670/J  Hopkins 
Olcy  Attn  DOC  CON  for  M  Sanpford 
Olcy  Attn  DOC  CON  for  R  Whittaker 
Olcy  Attn  DOC  CON  for  MS  364 
(Class  Reports  Lib) 

Olcy  *  Attn  DOC  CON  for  G  Spillman 
Olcy  Attn  DOC  CON  for  A  Davis 
Olcy  Attn  DOC  CON  for  R  Bridwell 

Lovelace  Biomedical  &  Environmental 
Research  Institute,  Inc. 

P  0  Box  5890 
Albuquerque,  NM  87115 

Olcy  Attn  R  Jones  (Unclas  only) 

Martin  Marietta  Corp. 

P  0  Box  5837 
Orlando,  FL  32855 

Olcy  Attn  G  Fotieo 

McDonnell  Douglas  Corp. 

5301  Bolsa  Avenue 
Huntington  Beach,  CA  92647 
Olcy  Attn  R  Halprin 

Merritt  Cases,  Inc. 

P  0  Box  1206 
Redlands,  CA  92373 

Olcy  Attn  J  Merritt 
Olcy  Attn  Library 


Lavrencu  Livermore  National  Lab. 

P  0  Box  808 
Livermore,  CA  94550 

Olcy  Attn  DOC  CON  fur  L-200  T  Butkevich 

Olcy  Attn  IXSC  CON  for  Tech  Infor  Dept.  Library 

Olcy  Attn  DOC  CON  for  L-205  J  Hearst  (Class  L-203) 

Olcy  Attn  DOC  CON  for  L-90  0  Norris  (Clsss  L-504) 

Olcy  Attn  DOC  CON  for  L-437  R  Schock 

Olcy  Attn  DOC  CON  fbr  L-7  J  Kahn 

Olcy  Attn  DOC  CON  for  L-96  L  Woodruff  (Class  L-94) 

Olcy  Attn  DOC  CON  fro  L-90  R  Dong 
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Meteorology  Research,  Inc. 

464  W  Woodbury  Road 
Altadena,  CA  91001 
Olcy  Attn  W  Green 

Nathan  M.  Newmark  Consult 
Eng  Services 
B106A  Civil  Eng  Bldg. 

University  of  Illinois 
Urbana,  IL  61801 

Olcy  Attn  N  Newmark 

Oak  Ridge  National  Lab. 

Nuclear  Division 
2-10  Lab  Records  Div 
P  0  Box  X 

Oak  Ridge,  TN  37830 

Olcy  Attn  Civil  Def  Res  Proj 
Olcy  Attn  DOC  CON  for  Central 
Research  Library 

Pacifica  Technology 
P  0  Box  148 
Del  Mar,  CA  92014 
Olcy  Attn  G  Kent 
Olcy  Attn  R  Bjork 

Physics  International  Co. 

2700  Merced  Street 
San  Leandro,  CA  94577 
Olcy  Attn  E  Moore 
Olcy  Attn  L  Behrmann 
Olcy  Attn  Technical  Library 
Olcy  Attn  F  Sauer 

R  &  D  Associates 
P  0  Box  9695 

Karina  Del  Rey,  CA  90291 

Olcy  Attn  Technical  Infer  Ctr 
Olcy  Attn  A  Utter 
Olcy  Attn  A  Kuhl 
Olcy  Attn  J  Carpenter 
Olcy  Attn  C  MacDonald 
Olcy  Attn  R  Port 
Olcy  Attn  J  Lewis 

Rand  Corp, 

1700  Main  Street 
Santa  Monica,  CA  90406 
Olcy  ATTN  C  How 

Sand  la  Laboratories 
Livermore  Laboratory 
P  0  Box  969 
Livermore,  CA  94450 

Olcy  Attn  DOC  CON  for  Library 
6  Security  Class  Division 


Sandia  National  Laboratories 
P  O  Box  5800 
Albuquerque,  NM  87185 

(Attn  Mail  Services  Section 
for  Intended  Recipient) 

Olcy  Attn  Mail  Ser  Sec  W  Roherty 
Olcy  Attn  Mail  Ser  Sec  3141 
Olcy  Attn  Mail  Ser  Sec  L  Vortman 
Olcy  Attn  Mail  Ser  Sec  A  Chaban 
Olcy  Attn  Mail  Ser  Sec  L  Hill 

Science  Applications,  Inc. 

P  0  Box  2351 
La  Jolla,  CA  92038 

Olcy  Attn  Technical  Library 

Science  Applications,  Inc. 

1257  Tasman  Drive 
Sunnyvale,  CA  94086 
Olcy  Attn  J  Dlshon 

Science  Applications,  Inc. 

2450  Washington  Avenue 
San  Leabdro,  CA  94577 
Olcy  Attn  D  Maxwell 
Olcy  Attn  D  Bernstein 

Science  Applications,  Inc. 

P  0  Box  1303 
McLean,  VA  22102 

Olcy  Attn  M  Knascl 
Olcy  Attn  B  Chambers  III 
Olcy  Attn  R  Slevers 
Olcy  Attn  J  Cockayne 

Southwest  Research  Institute 
P  O  Drawer  28510 
San  Antonio,  TX  78284 
Olcy  Attn  W  Baker 
Olcy  Attn  A  Uentel 

SRI  International 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025 

Olcy  Attn  0  Abrahaason 

Systems,  Science  A  Software,  Inc. 

P  O  Box  1620 
La  Jolla,  CA  92038 
Olcy  Attn  Library 
Olcy  Attn  D  Crlne 
Olcy  Attn  T  Rincy 
Olcy  Attn  X  Pyatt 
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Teledyne  Brown  Engineering 
Cummings  Research  Park 
Huntsville,  AL  35807 

Olcy  Attn  J  Ravenscraft 

Terra  Tek,  Inc. 

420  Wakara  Way 
Sclt  Lake  City,  UT  84108 
Olcy  Attn  Library 
Olcy  Attn  S  Green 
Olcy  Attn  A  Jones 

Tetra  Tech,  Inc. 

630  N  Rosemead  Blvd. 

Pasadena,  CA  91107 
Olcy  Attn  L  Hwang 
Olcy  Attn  Library 

TRW  Defense  &  Space  Sys  Group 
One  Space  Park 
Redondo  Beach,  CA  90278 
Olcy  Attn  I  Alber 
Olcy  Attn  Tech  Infer  Ctr 
Q2cy  Attn  N  Lipner 
Olcy  Attn  P  Bhutta 
Olcy  Attn  D  Baer 
Olcy  Attn  R  Plebuch 

TRW  Defense  &  Space  Sys  Group 
P  0  Box  1310 

San  Bernardino,  CA  92402 
Olcy  Attn  E  Wong 
Olcy  Attn  P  Dai 

Universal  Analytics,  Inc. 

7740  U  Manchester  Blvd 
Playa  Del  Rey,  CA  90291 
Olcy  Attn  E  Field 

Weidlinger  Assoc.,  Consulting  Eng 
110  E  S9th  Street 
New  York,  NY  10022 
Olcy  Attn  N  Baron 

Wetdlinger  Assoc.,  Consulting  Eng 
3000  Sand  Hill  Road 
Menlo  Park,  CA  94025 
Olcy  Attn  J  laenberg 

Westlnghouse  Electric  Corp. 

Marine  Division 
Hendy  Avenue 
Sunnyvale,  CA  94088 
Olcy  Attn  W  Volt 
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